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The Newman- Penrose formalism is used to derive the Teukolsky master equations 
controlHng massless scalar, neutrino, electromagnetic, gravitino, and gravitational 
field perturbations of the Kerr-de Sitter spacetime. Then the quantum entropy of 



a non-extreme Kerr-de Sitter black hole due to arbitrary spin fields is calculated 
by the improved thin-layer brick wall model. It is shown that the subleading order 
contribution to the entropy is dependent on the square of the spins of particles and 
that of the specific angular momentum of black holes as well as the cosmological 
constant. The logarithmic correction of the spins of particles to the entropy relies 
on the rotation of the black hole and the effect of the cosmological constant. 
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I. INTRODUCTION 



Ever since Bekenstein jl| and Hawking [2] discovered that black holes are thermodynamic 
objects endowed with a temperature proportional to its surface gravity and an entropy equal 



to one fourth o: 



its surface area, much effort has been devoted to investigating the statistical 



^ , quantum [4 1 , or dynamic origin of black hole entropy ^ . A proposal for the study of 
statistical origin of the black hole entropy is the brick wall method (BWM) suggested by 't 
Hooft 3|. In this model, the black hole entropy is identified with the statistical-mechanical 
entropy of a thermal gas of quantum field excitations outside the event horizon, which is 
composed of the leading order correction (i.e., the standard Bekenstein-Hawking entropy) and 
the logarithmic contributions to the black hole entropy. The subleading order corrections 
contain two parts, in general: the logarithmic term from the integral in the optical space 
and the logarithmic correction from the effective potential. This method, especially with 
the aid of the Newman-Penrose formalism ^ ^, has been successfully used in studies of the 
statistical- mechanical entropy of scalar fields for some static black holes 7, ^ and stationary 
axisymmetric black holes Q]. In the case of various static spherically symmetric black 
holes, it has also been applied to evaluate the entropy of spinor fields IJ, llJ| and electro- 
magnetic fields ^|. Recently the BWM has been used to calculate the entropy of Dirac field 
for the Kerr(-Newman) black hole However all these calculations did not consider the 
logarithmic contribution from the coupling of the spin of particles with the rotation of black 
holes. This kind of spin-rotation coupling effect appears [l^] in the Hawking thermal radi- 
ation spectrum of Dirac particles and photons in some non- stationary Kerr (-Newman) black 
holes 



18| . Though the subleading order correction had been included in Ref. but only 
the logarithmic contribution from the integral in the optical space had been considered there, 
the logarithmic term from the effective potential including the quadratic spin terms still had 
been neglected in those studies. The latter correction cannot be ignored in general, because 
it is of the same order as the former in the high frequency approximations. The presence of 
a logarithmic divergence in the entropy of the quantum scalar field has also been confirmed 
by other different approaches 3, 2^, 21 1, such as the conical singularity approach, the heat 
kernel expansion method, and the ^-function regularization technique, etc. As far as the ro- 
tating black hole case is concerned, the logarithmic contribution both from the integral in the 
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dragged optical space and from the effective potential had been considered by using the BWM 
only for the scalar field in recent research 
Recently much attention 



higher spin fields 



has been paid to the quantum entropies of black holes due to 
25 . 2^. Li, Shen and Gao, and Gao and Shen investigated 
the entropies of arbitrary spin fields in various spherically symmetric black holes, but did not 
consider the logarithmic contribution from the integral in the optical space nor the subleading 
order correction from the effective potential including the quadratic spin terms. Li studied 
the entropy of one component of massless fields with spin s = 1/2, 1, and 2 in the Reissner- 
Nordstrom black hole by the BWM and found that the logarithmic correction to the entropy 
depends on the linear term of the spins of the particles. Jing and Yan [25.] calculated the 
entropy up to subleading terms of massless fields of spin s > for the Kerr black hole and 
showed that the contribution of the spins to the logarithmic terms shall decrease the statistical- 
mechanical entropy of a Kerr black hole. Lopez- Ortega [2^ extended this analysis to the 
Rarita-Schwinger field case but pointed out that the entropy is increased by the logarithmic 
terms relating to the square of spins of particles. However the coefficients of the quadratic 
spin terms in the expressions of entropy presented in Refs. 3, 2^ are incompatible with each 
other and their validity is in doubt. Thus it deserves to take the logarithmic terms to the 
black hole entropy from the effective potential into account in details. How the spins of the 
quantum field changes the quantum entropy of a rotating black hole is an interesting question 
and needs to be further clarified. The knowledge of the entropy as a function of the spin of 
the field will be helpful to study the species dependence problem of the BWM on a rotating 
black hole. 

On the other hand, although the original BWM has contributed a great deal to the under- 
standing and calculation of the entropy of a black hole, there are some drawbacks in it such as 
the little mass approximation and taking the term including {L being the "infrared cutoff" ) 
as a contribution of the vacuum surrounding the black hole, etc. The model is constructed 
on the basis of thermal equilibrium at a large scale, so it cannot be applied to cases out of 
equilibrium, such as spacetime with two horizons, for example, a Schwarzschild-de Sitter black 
hole and Vaidya black hole {27, 2^. However one can improve this original BWM by taking 
only the entropy of a thin layer near the event horizon of a black hole into account, and utilize 
this improved thin-layer brick wall method to resolve some thermal nonequilibrium problems 
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that can hardly be treated by the original BWM 



27 



For a Kerr-de Sitter spacetime with non-degenerate horizons, it has a cosmological horizon 



and an outer black hole event horizon. As shown in Refs. [29|, |30|, there will appear two 
cases: (1) the general case that the temperatures of these two horizons are distinct. In this 
case, a non-extreme Kerr-de Sitter spacetime is a thermal nonequilibrium system. When the 
cosmological constant is very small, the two horizons will separate far away. Then each horizon, 
in principle, can be treated as an isolated thermodynamical system. Although the total system 
consisted of the two horizons is thermal nonequilibrium, the thin layer near the horizon can 
be taken as a local thermal equilibrium system. The quantum entropy of such a black hole 
can be calculated via this improved BWM which means that the entropy comes from a thin 
layer near the horizon. By arguments on thermodynamics of a composite system, the total 
entropy of a non-extreme Kerr-de Sitter spacetime is then taken the sum of the contribution 
from each horizon; and (2) a special case that the temperatures of the two distinct horizons are 
equal. The metric in this special case is called the lukewarm solution 29, 3^. The lukewarm 
solution achieves thermodynamic equilibrium, however it is unstable to changes in the mass 
of the black hole. In the lukewarm case, the notion of local thermal equilibrium may not be 
used, and one can still work with the original BWM to calculate the entropy of each horizon. 
It should be pointed out that the expression of the entropy in the lukewarm case is just a 
special case of that obtained in our general considerations. So we shall mainly consider the 
general case and leave the lukewarm case for a special discussion. 

In this paper, the entropies of non-extreme Kerr-de Sitter black holes due to higher spin 
fields are taken up for consideration on which the effect of the cosmological constant and that of 
the spins of particles are emphasized. The subject is important because the study of spacetimes 
which are asymptotically de Sitter has received a great deal of attention recently. The recent 
astrophysical observations of type la supernovae which indicate a positive cosmological 
constant 3^ and the recent dS/conformal field theory (CFT) correspondence js^ are two main 
motivations for studying the thermodynamics of black holes with a cosmological constant. It 
should be noted that a realistic black hole may be in an asymptotically nonflat space, thus it 
becomes important to investigate the effect of the cosmological constant on the entropies of 
these kinds of black holes. In Ref. js^, the renormalized black hole entropy for the massive 
scalar field in a Schwarzschild-anti de Sitter space has been considered recently by Winstanley 



via the original "brick wall" method. Although recent work [35[ has dealt with the entropy due 
to massless Dirac fermionic and scalar fields in the Newman-Unti- Tamburino (NUT) Kerr- 
Newman-de Sitter spacetime case, it puts emphasis on the improved brick wall model and only 
cares for the leading correction to black hole entropy, i.e., the standard Bekenstein-Hawking 
entropy. 

The purpose of this paper is to deduce expressions of the entropy of non-extreme Kerr-de 
Sitter black holes arising from arbitrary spin fields by using the improved brick-wall method 
and to investigate effects of the spins of particles and the cosmological constant on the statis- 
tical entropy. In this study, we carefully deal with the subleading order contribution to the 
entropy not only from the integral in the dragged optical space but also the logarithmic term 
from the effective potential including the quadratic spin terms, namely, the subleading correc- 
tions to the entropy arising from the coupling of the spins of the particles with the rotation 
of the black holes and the cosmological constant, regardless of a positive or negative one. 

The paper is organized as follows. Within the Newman-Penrose formalism P, 0] , we de- 
rive in Sec. II the master equations governing massless Klein-Gordon scalar, Weyl neutrino. 
Maxwell electromagnetic, Rarita-Schwinger gravitino, and linearized Einstein gravitational 
field perturbations of the Kerr-de Sitter space. Section III is devoted to deducing expressions 
of the statistical-mechanical entropy of the non-extreme Kerr-de Sitter black hole due to arbi- 
trary spin fields by using the thin-layer BWM. In Sec. IV, we give some arguments about the 
difference with previously published results. The last section summarizes our discussions. Ap- 
pendix A separates the Teukolsky master equations and presents the Teukolsky-Starobinsky 
identities in the Kerr-de Sitter geometry. In Appendix B, we calculate some integrals by the 
thin-layer BWM. 



II. PERTURBATIONS OF SPIN FIELDS IN THE KERR-DE SITTER SPACE 



The line element of the Kerr-de Sitter spacetime can be written in a coordinate system of 
Boyer-Lindquist type as 

HQ 



as = ^ [at — a sm ddip ) + 



adt — (r + a )dif 



l2 fdP de^\ , , 
Sf — + — ), (11 
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where 



A, = + a^) (^1 _ ^) _ 2Mr , A, = 1 + ^ cos^ 9, x = 1 + ^ , 



PP 



+ cos^ , p = r + ia cos 9 , p* = r — ia cos 6* . 



Here M is the mass of the black hole, a is its angular momentum per unit mass, and A = 3/1 



is the cosmological constant. To fit with the recent astrophysical observations [31| of type la 
supernovae, we require that the constant A is positive but very small j^^. This means that 
the radius of the de Sitter horizon is very large. The metric determinant, Ricci scalar, and 
the only nonvanishing Weyl scalar are, ^/—g = Ssin^/x^, TZ = AA = 12//^, \E'2 = —M/p*^, 
respectively. 

The metric has coordinate singularities at the roots of A,. = 0. For a range of parameters 
that satisfy the condition 1^/3 ^ > a^, there are four distinct real roots, three of which are 
positive and correspond (in decreasing order) to the cosmological, outer, and inner horizons of 
the black hole, respectively. The fourth root is negative and non-physical. A significant case 
is the lukewarm solution 2^ 3 which satisfies the condition = a^x^. In this special case 
the horizons are generally distinct whereas the surface gravity on both horizons is identical. 
In the allowed range of the physically meaningful parameters a and M, the admissible classes 
of the Kerr-de Sitter solutions are categorized completely by Booth and Mann js^ in details. 
In this paper, we shall assume that the horizons are non-degenerate and are mainly interested 
in the general case that the temperature of the cosmological horizon and that of the outer 
black hole event horizon are different from one another, and separate the lukewarm case for a 
special discussion. 

To derive a single master equation governing the perturbations of the Kerr-de Sitter space- 
time, we work it within the Newman- Penrose formalism P, |^ by choosing the null-tetrad 

vectors as 



1 

2S 



V2A~ep 



ixa sin 9, 0, 



sin 9 



(2) 
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and obtain the nonvanishing spin coefficients as follows: 

_ -^r- _ K _ -?av/A^sin6' 

where a prime denotes the partial differential with respect to its argument. 

Assuming that the azimuthal and time dependence of the perturbed fields will be of the 
form e*(™'^~'^*), we find that the directional derivatives are js^] 

6 = m^d, = ^Cl 6 = m'^d, = ^^Co, (4) 



where 



^" " 9r + A. + ""a. ' " 9^ + A. + ^1"°^^ 



and 



2A, 



Ki = ujir^ + a^) — ma , ^2 = 01^ sin 6 (5) 

sm 



with the relations 

Ki — sin 6 = tuS , i^^^ + i^r2 cot 9 = 2auj cos 9 . (6) 

Using the Newman-Penrose formalism P, 0] it can be shown that perturbation master 
equations in the Kerr-de Sitter geometry are separable for massless spin s = 0, 1/2, 1, 3/2, and 
2 fields 3^. The Teukolsky's master equations 39] controlling the perturbations of Kerr-de 
Sitter black hole for massless arbitrary spin fields (s = 1/2, 1, 3/2, and 2 for Weyl neutrino, 
source-free Maxwell electromagnetic, Rarita-Schwinger gravitino, and the linearized Einstein 
gravitational fields, respectively) read ji^ 

I [D - (2s - l)e + e* - 2sp - p*] (A - 2s7 + /i) 
-[S- (2s - 1)P -a* - 2sT + n*] {6 - 2sa + vr) 

-(s-l)(2s-l)^2|$s = 0, (7) 
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for spin weight s = 1/2, 1, 3/2, 2 and 



I [A + (2s - 1)7 - 7* + 2s/i + /i*] {D + 2se - p) 
-[6 + (2s - l)a + + 2s7r - r*] {6 + 2s(3 - r) 
-(s-l)(2s-l)^2}$-s = 0. 



(8) 



for spin weight s = —1/2,— 1,-3/2,— 2 (s is the spin of the perturbed fields). It is clear that 
Eqs. ((Zj) and (jH)) are also valid when s = 0, they coincide with the massless minimally coupled 
scalar field equation 

(□ + 7^/6)<l> = 0, (9) 

with $ = $0 = '^'-o- 

All the above equations are separable by using the Newman- Penrose formalism, and can 
be written as (ignoring the factor e*'-'"'''"'^*-') 



{ArV^Vl + ^eCU^eCs) + 2{2s - 1) - 



= 0, (10) 
(p*2^$_,) =0. (11) 



It can be directly shown that Eqs. (fTUj) and (ITIT) a re also satisfied by the scalar Debye potentials 
<Ps = ^s/p*"^"" and = p*^''$_s, which obey j4o| 

I - (2s - l)e + e* - p*] [A - 2s7 - (2s - l)p] 
-[5- (2s - l)p -a* + n*]\6- 2sa - (2s - l)7r] 

-(s-l)(2s- 1)^2)0^ = 0, (12) 
I [A + (2s - 1)7 - 7* + p*] [D + 2se + (2s - l)p] 
-[6 + (2s - l)a + 13* - T*] [6 + 2s(3 + (2s - l)r] 

-(s-l)(2s- 1)^2)0-^ = 0. (13) 

In Ref. js^, the authors showed that Eqs. (|T0|) and (jllll can be further separated and 
transformed into Heun equations. In the case of a Kerr black hole, they degenerate to the 



generalized spheroidal wave equation 







42| . a confluent form of Heun equation [42]. Exact 



solutions to these equations, and i nteg ral equations as well as other related applications can 



be found in Refs. 
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42 



44l |. To the end of this paper, we do not take these into 



account, but just present in Appendix A the radial Teukolsky-Starobinsky identities 



40 



where the coefficient C2 corrects the previously pubhshed results 
expressions of the above equations (fTUj) and (fTT|) 



From their obvious 



Si " arV " dr 



sm 9o9\ oO 



1 



A' 



XK2- s[^ + Ae cote 







(14) 



and 



i|A^^fA^-^) + 



X^Kl + isxK^K s 



1 a 



A. 



- -A" + 

2 " sin^a^ 



Ag sin 



1 

'A^ 



A' 



xK2 + s[^ + Aecoie 



Aisxojp 



4g^ + 2, 

/2 ' 



0, 



(15) 



one can easily find that they are dual by interchanging s — —s. Thus one only needs to 
consider the case of positive spin state p = s, and obtain the results for the negative spin 
state p = — s by substituting s ^ —s. Eqs. dH)) and (fT3j) can be combined into the form of 



Teukolsky's master equation 



A, 92 ^ (1 + s)a: d 



S 9r2 ' S dr 



A^dP_ 

S 



2^,2 



A'q + Ae cot 6' d ^ uj^x 



89 



-2 + a^V sin^ 9 



A, 



S 

iA' 



Ar 



(r^ + a^) + 2ir 



1 \ m?x^ ( 
A^J ^ S Ia^ 
iaAL 



2A, 

s^A, / A' 



2smx 



2A. 



Ae sin^ 9 
sine V2A, 



2sc^ 
S 

cote 



a sm 



/ A' 
9{ — ^ - cote 



V2A 



^A" 

2S ' 



4s2 + 2 



2Afl 



+ cote = 0, 



[s = 0, ±1/2, ±1, ±3/2, ±2). 



(16) 



In the above master equation (jl6p . one can identify the term proportional to (4s^ ± 2)//^ = 
(2s^ ± 1)7^/6 with the conformally coupled term for arbitrary spin fields in virtue of the 
nonvanishing Ricci scalar in the Kerr-de Sitter spacetime. In the next section we will utilize 
this equation to obtain the density of states by using the Wentzel-Kramers-Brillouin (WKB) 
approximation scheme. 



III. ENTROPY OF KERR-DE SITTER BLACK HOLES DUE TO SPIN FIELDS 



Now we calculate the entropy due to arbitrary spin fields for the non-extreme Kerr-de 
Sitter black hole by the thin-layer BWM. First we try to seek the total number of modes 
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with energy less than uj. In order to do this, we make use of the WKB approximation and 
substitute ~ e**^*^'''^^ into the above Teukolsky's master equation (fTBj) . then we obtain 



2 2 

m X 



1 



+ 



Ae sin^ 9 A, 
s^A, / A' 



cot 6* 



cote 



A^ 
2A. 



A, Ae 
2smx f A'q 



Ssin^ V2Afi 



cot 6* 



-—A" 

22 *^ 



0, 



:i7) 



S V2Ae / /2 

in which, kr = G^r and kg = are the "wave numbers." In terms of the covariant metric 
components g^^, Eq. (fTTj) can be rewritten as 



k 



ka 



grr gee 



g^^uj"^ + 2gt^muj + gum^ 
V 



+ 2{ujB + mG) +H, = 0, 



[18) 



where 

9w 



9tt 



E _ S 

a;' ^''-Ae 
(r2 + a2)2A« - A,a2 sin^^ 



Aga"^ sin^ 6 - Ar 



9t^ 



Ar - (r^ + a2)Ae 



.Xasing. 

S V 
s^A, / A' 



sin^ ^ , V = gttg^^ - gj^ = 
A^ \ ^ / a;, 



2A, 



2 45^ + 2 



Ssin^ V2A, 
-A'' . 



A,. Ae sin^ 6 

X' 

+ cot e) , 



a sin 6' 



(19) 



E V2A0 ■ --^v ■ /2 2E" 
To obtain the density of states, let us suppose that the quantum field is rotating with an 
angular velocity Qh in a thin-layer very near the horizon of the non-extreme Kerr-de Sitter 
black hole. Substituting S = uj — mQh into Eq. (jl8|) . we reduce it into the form 



9rr 9ee 



+2 



SB + m{nhB + C) 



+ Hs = 0. 



and then rewrite it as 



kn 



+ — + 

9rr gee 



-9tt 
-V 



m + mo 



-9tt V ^ 



(20) 



(21) 



where 



9tt = gtt + ^gt^Qh + 9w^h 5 

{gt^ + g^^nh)S + {^hB + G)V 



mo 
sW 



9tt + '^gtifi^h + 9<f<f^l 



{9tt + 9t<p^h)B - {gt^ + g^^Vth)G , 
Vs = H,- {guB^ - 2gt^BG + g^^G^) = P 



sA'; 

2S 



(22) 
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Qtt is the temporal component of the metric of the dragged optical space j47|, 4^, is the 
angular velocity caused by the rotation of the black hole and can be called the "spin potential," 



while Vs is the effective potential 



ll\ which is equal to /i^ + ^71 = /i^ + 12^ in the case 



of a massive scalar field with an arbitrary conformally coupling = 1/6 for the minimally 
coupling), where /i is the mass of the field. 

For a given energy u = S + mQh and a given azimuthal angular momentum m, Eq. 1)211) 
represents an ellipsoid of three-dimensional momentum space, a compact surface, spanned 
by kr, kg, and m, which is a subspace of a six-dimensional phase space, supposed that the 
following conditions could be satisfied: 

{£ + sW)'' + gttVs > . (23) 
ds the number of modes with S is 



t/rr > , gee>0, -gtt > , -V > 0, 



Therefore in this case, for the positive spin state of spin fie 
equal to the number of states in this classical phase space 

1 



T{£,s) 



(27r)3 
1 



3n 



drdOdif I dkrdkgdm 



de 



dr- 



-9 



(S + sWy + guVs 



1 3/2 



(24) 



'r-h+e [~9tt) 

Here we impose the improved thin-layer BWM boundary conditions in the integral with 
respect to the radial coordinate r, that is, we assume that the quantum field is equal to zero 
for r < rh + e and r > Vh + Ne, with rh ^ £■ The ultraviolet cutoff e is a small distance 
from the horizon to the inner brick wall, the cutoff parameter iV is a sufficient big integer 
to remove the infrared divergence. In other words, the infrared cutoff L in the original BWM 
is now replaced by Vh + Ne in the upper limit of the radial integral. This reflects a signiflcant 
difference from the original BWM, and such improved BWM is called the thin-layer model 
[2^, it can overcome some defects in the original BWM. 

In the original brick-wall model, it is supposed that the black hole is in thermal equilibrium 
with the external fleld in a large spatial region. This method cannot be applied to a nonequi- 
librium system such as a system of spacetime with two horizons, for example, Schwarzschild-de 
Sitter spacetime and Vaidya spacetime 3, 2^- In such spacetimes, two problems arise: (a) 
the thermal equilibrium between the external fleld and the hole is unstable, so the thermal 
equilibrium on a large scale basis for the brick wall does not exist; and (b) since the two hori- 
zons have different temperatures, there exists no global thermal equilibrium over the entire 
spacetime, and statistical physics laws are invalid there. 
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In the thin-layer method, the entropy of the black hole is mainly attributed to the degrees 
of the freedom of the field in the thin-layer {vh + e < r < + Ne) covering the surface of the 
horizon. The Bekenstein-Hawking entropy should be associated with the fields in this small 
region near the horizon, where the local thermal equilibrium exists and statistical physics laws 
are still valid. To guarantee that the notion of local thermal equilibrium can work very well, 
here one supposes that the physical quantities of the thermodynamic properties of the exciting 
field outside the hole vary slightly in the vicinity of the horizon. On the one hand, the length 
of the thin-layer region near the horizon must be small enough on a macroscopic scale so that 
the physical quantities in the region can approximately be treated as some constants in the 
vicinity, and approximate equilibrium in the small region is achieved. On the other hand, 
the region must be large enough on a microscopic scale so that the statistical mechanics are 
valid for the fields in the near horizon region, and the thermodynamic variables can be defined 
through a partition function. In order for local equilibrium to be maintained, it is necessary 
that the hole's radiation is slight enough that the fluctuation of thermodynamic properties of 
fields can be treated as a small quantity. In most cases, this condition can be satisfied except 
for those ones at the Planckian scale. 

In the spacetime that has two horizons with two different temperatures, there exists no 
global thermal equilibrium in the entire spacetime, however approximate thermal equilibrium 
exists in the two layers near each horizon. Thus the global thermal equilibrium is not needed, 
the validity of local thermal equilibrium is crucial to the discussion j28|. In such a thin-layer 
BWM, the total entropy is mainly attributed to the two thin-layers near the two horizons, 
namely, it is a linear sum of the area of each horizon. 

In the case of a non-extreme Kerr-de Sitter black hole with non-degenerate horizons, situa- 
tion is more involved. The issue of local versus global thermal equilibrium is a rather delicate 
one. As mentioned before, there are two cases that can happen: (1) In the general case, 
the temperatures of the cosmological and black hole event horizons are distinct; and (2) in 
the lukewarm case, both temperatures are equal. In the general, non-extreme Kerr-de Sit- 
ter background, there is no time-like Killing vector which is well-defined in the whole region 
surrounded by the two horizons. Hence the system in the general case cannot be in thermal 
equilibrium globally, namely, there exists no thermal equilibrium in the entire spacetime since 
the two horizons have different temperatures. In addition, there cannot exist a global ther- 
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mal equilibrium between the external field and the hole in a large spatial region. Thus in 
this general case, one must work with the thin-layer method to calculate the entropy of the 
Kerr-de Sitter black hole. In the lukewarm case (M^ = a^x^), however, the cosmological and 
black hole event horizons are in thermal equilibrium (but unstable to changes in the mass of 
the black hole 2.9] ). In this special case, one may not introduce the notion of local thermal 
equilibrium and still adopt the original BWM to calculate the entropy of each horizon. The 
main interest of this paper is in the general case, and we will turn to the lukewarm case for a 
special consideration. In both cases, the total entropy of a quantum field in the non-extreme 
Kerr-de Sitter spacetime is summed up from the entropies corresponding to the maximal and 
minimal spin-weight components. 

Summing over the positive and negative spin states p = ±s, we get the total states number 



bir 



Th+Ne 



dr- 



-9tt) 



1 3/2 



dr- 



-9 



9_ 
SttL 



Tu+e \-9tt 



(S + sWf + guVs 



i£-sWf + gttV. 



1 3/2 



4S 



-gttW 



(25) 



In the above, we have expanded Eq. ()25p in the high frequency approximation and introduced 
an appropriate degeneracy Qs for each species of particles (It is well-known that Qs = 2s + 1 
in the non-relativistic quantum statistics case). Here Qs = 1 for scalar field {s = 0), Qs = 2 
for Weyl neutrino {s = 1/2), Maxwell electromagnetic {s = 1), Rarita-Schwinger gravitino 
(s = 3/2) and linearized Einstein gravitational {s = 2) fields, and Qs = 4: for massless Dirac 
field {s = 1/2), respectively. The following four integrals in terms of the thin-layer BWM are 
calculated in Appendix B 



h 
h 



de 



r-h+Ne 



dr- 



-9 



9tt 
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Apparently the integral Ji is related to the volume of the dragged optical space 



(26) 
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It is known that a "physical space" must be dragged by the gravitational field with an 
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azimuthal angular velocity ilh in the stationary rotating spacetime. Since we have supposed 
that a quantum field in a thin-layer very near the horizon is in local thermal equilibrium 
with the non-extreme Kerr-de Sitter black hole at the asymptotic temperature measured 
by an observer located at the spatial infinity, it is appropriate to assume that the quantum 
field is rotating with angular velocity Qh in this thin-layer also. Note that here we use the 
asymptotic quantities rather than those measured by a local (corotating) observer in the 
thin-layer. In fact the equivalence principle implies that a thin-layer system in local thermal 



e,un,bH„,„ has a .oca. To.„,a. inve.se temperature give, by ^ p^., wit.> p being 

the asymptotic inverse temperature [42|. The temperature measured by this local observer 
is a local temperature Piocai, correspondingly he also measures a local (blue-shifted) energy 
i^iocai = ^/ \/ —9tti where the energy uj measured by the asymptotic observer is associated 



with the coordinate time t. The quantity (iuj = Piocai^^ local is an invariant due to the first 
thermodynamic law, and so is P{uj — mQh), etc. As the equivalence between using local 



quantities and using asymptotic ones has already been proven in Ref. j50|, one can directly 
adopt the asymptotic quantities, needless to make further conversions. For such a local (quasi- 
) equilibrium ensemble of states of spin fields, the free energy can be expressed as follows: 
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g/3(<^-mnh) _ ( — l)2s 

where we have made a variable substitution £ = uj — mVth- The formulas for the density of 
states g{(jj, m) = dT{u!, m)/duj and r(£^) = J V{£ + mVL^, m)dm had also been used. Inserting 
Eq. (j25|) into Eq. (j27|) and then carrying out the integral over we get the expression for 
the free energy 



F = 

6tt 



X ■ 



de 



Th+Ne 



dr- 



-9tt) 



d£ 



o(3£ _ (_ 



9s_ 

371 



{£ + sWf + guV, 
d£ 



3/2 



+ 



{£-sWf + guV-, 



3/2 



-9s 



eP£ _ (_i)2s 

\2s 



h£'' + 3{l2 + s^h)£-3s^h 
3 + (-p2s 



477/^2 



(28) 



15 



where ({n) = ^ 1//;;" is the Riemann zeta function, ({4) = 7r^/90, ({2) = 7r^/6, etc. 

k=l 

We are now ready to obtain the entropy of the non-extreme Kerr-de Sitter black hole due 



to arbitrary spin fields from the standard formula S = P'^{dF/dP), 

15 + 
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9s_ 
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-C(l) l + (-ir U^Z 



(29) 



Next, we are in a position to consider the four integrals Ji ~ 14. It is easy to obtain the 
Hawking temperature and the area corresponding to the horizon of the non-extreme Kerr-de 
Sitter black hole 



"27r 
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a: 



Ah = 4n{rl + a'')/x. 



(30) 



"^^Xirl + a2) x^^h ' 

By means of the thin-layer BWM, we take the angular velocity of a quantum field near the 
horizon of the non-extreme Kerr-de Sitter black hole as Qh = O'/i'^h + '^^)? find that only 
the integrals Ji and I2 contribute to the leading and subleading terms in the entropy, while 
the integrals I3 and I4 can be ignored as usual because the integral I3 can be attributed to the 
contribution of the vacuum surrounding the hole due to Is ~ 0{e) and /4 vanishes at least up 
to the order 0{e). The final expressions for these integrals are presented in Appendix B where 
the ultraviolet cutoff e is replaced by the proper distance t] from the horizon to the inner brick 
wall r] = j^^^^ yfg^dr 2{ellh/ ^'^^Y^"^ ■ In order to be comparable with the results already 
appeared in the literature and to simplify the expression, we have set the new ultraviolet cutoff 
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25 



e and infrared cutoff A by r/^ = 2e^/15 and = h?/e^ as did in Refs. 
utilization of the results given by Eq. ()B7|) . we obtain the statistical-mechanical entropy 
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Assuming that the field is in the Hartle-Hawking vacuum state and taking /3 = we get 
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that the entropy is given by 
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(32) 



Equation ()31|) or ()32|) shows that the entropy of a non-extreme Kerr-de Sitter black hole due 
to arbitrary spin fields consists of two parts, the leading order contribution and the subleading 
order corrections, or equivalently the contribution from the integral Ji in the dragged optical 
space and the logarithmic term I2 from the effective potential including the quadratic spin 
terms. The logarithmic corrections consist of the one from the integral Ji and that from the 
integral I2, both of them are of the same order, therefore the latter cannot be thrown away 
any more. It should be noted that the coefficients of the logarithmic divergence is universal 
(invariant under a change in the value of the cutoff, or even under a change in the regulator 
scheme). The expression of Eq. fl32|) may settle down the species dependence problem of the 
brick wall entropy on a rotating black hole, it covers many previously obtained results. 

To see what role the conformally coupling plays in the calculation, we also present the 
expression of the entropy for a massive scalar field in the non-extreme Kerr-Newman-de Sitter 
black hole 
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arh Vft/ 



e 



(33) 



with an arbitrary conformally coupling constant ^. It is not difficult to find that the confor- 
mally coupling will contribute a logarithmic correction to the entropy. By setting = 1 in 
Eqs. (1221) and (j^Sl), we recover the standard Bekenstein-Hawking entropy. 



IV. DISCUSSIONS 



The final result that we obtain for the entropies of arbitrary spin fields in the non-extreme 
Kerr-de Sitter space deserves some remarks. 
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(a) The entropies given above have summed up the contribution from the maximal and 
minimal spin- weight states of a quantum field. Under the condition that satisfies /^/3 ^ 

> a^, the cosmological horizon separates far away from the outer black hole horizon, the 
total entropy of the non-extreme Kerr-de Sitter black hole can take a linear sum of the area 
of the two horizons. 

The calculations here are valid both for the outer black hole event horizon case and for the 
cosmological horizon case because the two horizons are in an equal position. We think it is 
also valid for the black hole event horizon of the Kerr-anti de Sitter space by changing the 
sign of the cosmological constant A. 

(b) The entropies depend not only on the spins of the particles but also on the cosmological 
constant except different spin fields obey different statistics. They rely on the quadratic terms 
of and —l/t^ as well as a^. 

(c) Both the contribution of the spins and that of the cosmological constant to the entropies 
are in subleading order. The spins have a tendency to increase the entropies, but the effect of 
a positive cosmological constant tends to decrease them. 

(d) The logarithmic term from the spins of the particles not only depends on the spin- 
rotation coupling effect but also on the coupling between the spins of particles and the cos- 
mological constant. Figure 1 shows that how the coefficient of the square term of the spins 
for the logarithmic correction to the entropies is affected by the specific angular momentum 
of the hole and the cosmological constant. 

(e) Three special cases may be very interesting: 

Case I: In the Schwarzschild-de Sitter black hole case, the specific angular momentum 
vanishes and the black hole is non-rotating. The improved thin-layer model still works in this 
case, then the entropy due to arbitrary spin fields is 



where Ah — 47rr^. 

It is not difficult to find that the entropy still depends on the square of the spins of particles 
unless in the Schwarzschild and Reissner- Nordstrom black hole cases. The existence of a 
positive cosmological constant will decrease the entropy whereas the effect of a negative one 
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4 127r/2 




(34) 
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FIG. 1: The coefficient of the quadratic spin term for the logarithmic correction to the entropies is 
dependent on the specific angular momentum of the hole and the cosmological constant. The entropy 
increases when |a/r^| becomes large and cuts down when r^/l increases. 



tends to increase it, vice versa. This result is in accordance with the entropy of a massive 
scalar field for the Schwarzschild-anti de Sitter black hole by using the original BWM in Ref. 
3^ where a negative cosmological constant is assumed and if the square mass of the scalar 
field /i^ is replaced by 

Case II: In the Kerr black hole case, the cosmological constant A goes to zero (or in the 
/ oo limit), the improved thin-layer model degenerates to the original BWM by simply 
setting the infrared cutoff L = Ne without any change in the ultraviolet cutoff e. By the 
equality = A^/e^ and the relation e ^ r]"^ ^ e^, the infrared cutoff L A^e^ ~ A^ is 
consistent with the definition A^ = Le'^/e as given in Refs. 



201 3, Q- (There is a minor 



difference in the definition of the infrared cutoff L because A^ can be sufficient large, here it is 

n 

a shift by rh from the original infrared cutoff introduced by 't Hooft [7[). It should be noted 
that the system in the brick wall region is, in fact, still in a local thermal equilibrium, and the 
thin-layer model still works very well. Taking these into account, the expression of the brick 
wall entropy in the Kerr black hole case reduces to 

15 + (-1)2^ 
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where Ah = 4:iT{rl + a^). 

This resuh is very similar to that presented in Refs. 
coincides with the result given by Ref. j^]. In the cases of s = 1/2 (Dirac field), s = 1 



25 



26l. In the scalar field case, it 



(Maxwell field), s = 3/2 (Rarita-Schwinger field), and s = 2 (linearized gravitational field), 
it has the same coefficients as that in Ref. I25I. But these coefficient terms have a 
different tendency because they differ by the terms in the square bracket. In the Rarita- 



Schwinger field (s = 3/2 
as pointed out by Ref. 



case, our result manifests that the spins will increase the entropies 
2^. Again they differ by the coefficients of the quadratic term s^. 



The expression P = s^a? cos^ 6 [A 



S] /T? used in Refs. |2^ 



2611 but written in 



our notations, is probably the result that those authors had missed some terms in the process 
of their approximation. We point out that it should read P = As^a? cos^ ^ [A — (r^ + a^)] /S^, 
and these coefficients presented in Refs. j^, 2^ are incorrect in the Kerr black hole case. The 
difference of the coefficients of the quadratic term (the pure spin-dependent term) between 



ours and those in Refs. 



25 



26| is compared in Fig. 2. 




FIG. 2: Compare the coefficient of the pure spin-dependent term to tlie logarithmic correction. The 
curve denoted by J&Y ^| demonstrates that the entropy is decreased by the contribution from the 



term. Both the curve denoted by ALO 



26l | and that of ours denoted by W&Y illuminate that the 



spin-rotation coupling effect will increase the entropy of a Kerr black hole. 



Q, Q, 



Case III: In the lukewarm case 12a Ian, the non-extreme Kerr-de Sitter spacetime is 
characterized by the condition = a?"x^. In the absence of rotation, the lukewarm solution 
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reduces to de Sitter spacetime. When the cosmological constant vanishes, it degenerates to 
the extreme Kerr black hole. Obviously the lukewarm solution is a special case of what we 
have discussed in the above. In this case, the cosmological and outer black hole horizons are 
distinct, but the temperatures on both horizon are equal, so the cosmological and black hole 
event horizons are in thermal equilibrium. Just as the Kerr black hole case, here one need not 
assume local thermal equilibrium and he can still work with the original BWM to calculate 
the entropy of each horizon. Of course, the cutoff L cannot be so large as to close to the 
cosmological horizon. The final result of the entropy is still expressed by Eq. fj32|) with the 
special restriction (namely, = a^x^)on the allowed range of parameter of the lukewarm 
solution. It should be noted that the same result can also be arrived at by using the thin-layer 
method. This issue can be easily understood as the brick wall entropy is mainly attributed to 
the degrees of freedom of the field in the near horizon region. As such we think the thin-layer 
method is more universal than the original BWM. 

(f) We do not consider here the contribution from the super-radiant modes. As is well 
known, for rotating black holes in an asymptotically flat space, classical super-radiance effects 
occur only for bosonic but not fermionic fields however the quantum analog of super- 
radiance does occur for both bosonic and fermionic fields l51|. It has been pointed out that 

n 

the super-radiant modes have a contribution to the entropy of a rotating black hole 1521] . In the 
Kerr black hole case, the effect of the super-radiant modes is to halve the leading contribution 
to the entropy of scalar and Dirac fields jl^ . \i\ . Their affect on the subleading term to the 
entropy for the non-extreme Kerr-de Sitter black hole is an interesting problem and deserves 
to be discussed elsewhere. 

(g) We do not address here the renormalization of the divergence in the entropy of the non- 
extreme Kerr-de Sitter black hole. A generally accepted belief is that the matter field contribu- 
tions to the entropy can be interpreted as one-loop corrections to classical Bekenstein-Hawking 
entropy. The divergence of the brick wall entropy has the same origin as the divergence of 
the one-loop effective action in quantum field theory in curved space. The leading quadratic 
divergence of the entropy calculated by the brick wall model can be absorbed by the renor- 
malization of the Newton gravitational coupling constant in the one-loop effective action for 
matter and gravitation fields 53]. When gravity is described by a higher curvature effective 
action, the standard Bekenstein-Hawking result is only the leading contribution in the entropy. 
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and there are still additional corrections from the higher curvature interaction js^. On the 
other hand, quantum corrections in curved s pac e are known to result in higher order curvature 
contributions to the Einstein-Hilbert action [55| . The subleading order logarithmic divergence 
in the entropy requires the introduction in the gravitational action of term quadratic in the 



curvature to be renormalized 



56l | (see also 
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57 



). The standard renormalization of 



Newton gravitational coupling constant, the cosmological constant and other coefficients by 
the curvature squared (i?^) terms in the one-loop effective gravitational action 55| can re- 
move all the divergent (quadratic and logarithmic) terms in the entropy 5^ , so the remaining 
quantity is finite. 

In the case of non-extreme Kerr-de Sitter geometry, it is clear that the divergence coming 
from the integral in the dragged optical space can be completely removed by the above pro- 
cedure, but it is unclear to us whether the logarithmic divergence arising from the effective 
potential including the quadratic spin terms needs further introduction of other interaction 
terms to be renormalized because this divergence explicitly contains quantum corrections to 
the brick wall entropy from the spin-rotation coupling interaction. Though the renormaliza- 
tion of the divergence in the brick wall entropy of the non-extreme Kerr-de Sitter black hole 
has not been discussed here, we hope the expression of the entropy obtained here can shed 
light on this subject because there is still little work on dealing with this problem in a rotating 
black hole spacetime. 



V. CONCLUSIONS 

In summary, we have studied the statistical entropies of the non-extreme Kerr-de Sitter 
black holes due to arbitrary spin fields, especially the subleading corrections to the black 
hole entropy arising from the coupling of the spins of particles with the rotation of the black 
holes. First, the null-tetrad in the Newman-Penrose formalism is introduced to decouple the 
Teukolsky master equations governing massless scalar, neutrino, electromagnetic, gravitino, 
and gravitational field perturbations of the Kerr-de Sitter space. Then, starting from the 
Teukolsky master equations we seek the total number of the modes of the fields by taking the 
WKB approximation. Last, the free energy and the quantum entropy of a non-extreme Kerr- 
de Sitter black hole due to arbitrary spin fields are calculated by the improved thin-layer brick 



22 



wall model. It is shown that the subleading order contribution to the entropy is dependent 
on the square of three quantities, the spins of particles, the specific angular momentum of 
black holes, and the cosmological constant. The contribution of the spins of particles to the 
logarithmic terms of the entropy depends on the spin-rotation coupling effect and the effect 
of the cosmological constant. It should be noted that the final result also holds true in the 
lukewarm case where local thermal equilibrium need not be assumed. 

In particular, we have carefully investigated the effect of a positive cosmological constant 
on the black hole entropy and shown that the correction from the effective potential is of the 
same order as that from the integral in the dragged optical space, and both of them cannot be 
discarded. However it should be noted that the correction from the "spin potential" W has 
been attributed to that of the vacuum surrounding the black holes and it has been neglected 
here. Possible new quantum effects related to the "sW" term in a Kerr-de Sitter black hole 
geometry may be a very interesting thing and deserve to be further investigated. It is also 
needed to extend this analysis to the Kerr-Newman-de Sitter black hole case. 
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APPENDIX A. TEUKOLSKY-STAROBINSKY IDENTITIES 



The Teukolsky's master equations (fTUI) and governing the perturbation of Kerr-de 
Sitter space with massless fields can be separated into the angular parts 
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and the radial parts 
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Here we only consider the radial equations. It is clear that {Af.Rs) and R-s are proportional 



to complex conjugate functions. The exact re^ 



famous Teukolsky-Starobinsky identities 
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ationship between these functions are called the 
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with the coefficients 
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where a"^ = — ma/uj. The coefficient IC2P corrects the previous results |46'], while the 
appearance of IC3/2P is the first time, to our knowledge. 



APPENDIX B. INTEGRALS IN TERMS OF THIN-LAYER BRICK WALL MODEL 



Distinguished from the original BWM, the thin-layer BWM suggests that the entropy of a 
black hole with two horizons mainly comes from a very thin layer in the vicinity of the horizon 
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where exists a local thermal equilibrium. Just as the original BWM, it also impose a small 
ultraviolet cutoff e such that 



^(a;) = for r < r/j + £ , 



(Bl) 



where denotes one coordinate of the two horizons of the non-extreme Kcrr-dc Sitter black 
hole. In this paper, it represents the outer black hole event horizon or the cosmological horizon, 
satisfying the horizon equation 

„2 , 



A,, = (r^ + a2)(l-^)-2Mr, = 0. 
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To remove the infrared divergence, it however introduces another cutoff parameter — an 
arbitrary big integer N such that 



= for r>rh + Ne. 
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Suppose that the quantum field is rotating with the angular velocity Vlh = a/irl + a^) in 
the thin layer near the horizon of the Kerr-de Sitter black hole, we may expand A,, close to 
the event horizon rh as 
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and then expand three quantities gu, P, and W in terms of the surface gravity Kh 
KJiMrl + a')) = 27rA'J(x'An) as follows: 
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where E^ = + cos^ ^. 
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Expanding the integrands in the four integrals Ji ~ I4 defined in Eq. 
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and carrying out the integrals with respect to 9 and r, we finally arrive at 
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where = 'iTi{r'j^ + a^) / x is the horizon area. In the last step, we have replaced the ultraviolet 
cutoff e = rj'^A'^ /{AHh) by the proper distance rj from the horizon to the inner brick wall 
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fr ^fg^rdr ^ 2{e^h/KhY'^- To be comparable with Refs. 



3, 



25 



26|, the new 



infrared cutoff A and ultraviolet cutoff e in the above equation ()B7|) are defined by 
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For large A^ and small a, Eq. ()B7|) implies 772 = 2e2/15, which is used in the context. 
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